ifcWWP35«#!^^  ir%W^«!i^^ 


AD/A-006  391 


FORMULATIONS  OF  THE  EQUATIONS  OF 
MOTION  OF  AN  ELASTIC  AIRCRAFT  FOR 
STABILITY  AND  CONTROL  AND  FLIGHT 
CONTROL  APPLICATIONS 

Robert  C,  Schwanz 

Air  Force  Flight  Dynamics  Laboratory 
Wright-P.atterson  Air  Force  Base,  Ohio 

August  1972 


J 


NifilS 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


AD  A 0 06391 


FORMULATIONS  OF  THE  EQUATIONS  OF  MOTION 

OF  AN  ELASTIC  AIRCRAFT  /? 

FOR  STABILITY  .AND  CONTROL  . /'* 

AND  FLIGHT  CONTROL 
APPLICATIONS  . 


/ 


Robert  C.  Schwanz 
Control  Criteria  Branch 
Flight  Control  Division 


August  1972 


Rep'odu<*ed  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


u$  Department  of  Commorco 
Springfioid,  VA.  2215t 


t 


Air  Force  Flight  Dynamics  Laboratory 
Air  Force  Systems  Command 
Wright-Patterson  Air  Force  Base,  Ohio 

riSihTKl 
1 Apprt 


rjUliOH  S'L AT L 


>vt*d  lor  public 


% 


AD  AO  06391 


Abstract 


This  ^Technical  Memorandum  presents  a derivation  of  the  equations  of 
motion  for  large  and  small  disturbance  perturbations  from  a reference 
state  of  motion.  The  small  perturbation  equations  of  motion  are  then 
{ generalized  to  include  the  effects  of  atmospheric  turbulence  and  gusts 
on  the  controls-free  elastic  aircraft.  The  resulting  equations  are 
termed  the  EXACT  formulation.  These  EXACT  equations  may  be  used  to 
evaluate  the  stability  and  performance  of  integrated  flight  control  sys- 
tems of  the  Control  Configured  Vehicle  (CCV)  type.  The  equations  are 
unique  in  that  they  describe  a highly  damped  system  using  motion  coordi- 
nates referenced  to  a body  fixed,  "mean",  non- inertial  axis. 

The  EXACT  formulation  of  the  equations  of  motion  are  difficult  to 
solve  numerically.  This  difficulty  can  be  overcome  in  some  flight  control 
analyses  by  using  a simplified  formulation  such  as: 

• QUASI  STATIC 
•MODAL  SUBSTITUTION 
•RESIDUAL  STIFFNESS 
•RESIDUAL  FLEXIBILITY 
•MODAL  TRUNCATION 

The  range  of  the  applicability  of  each  of  the  specialized  fonaulaticnc 
is  limited  by  the  assumptions  required  to  reduce  the  EXACT  formulation  to 
that  specialized  formulation.  A discussion  of.  the  assumptions  is  pre- 
sented to  guide  the  application  of  each  formulation  to  military  and  com- 
mercial aircraft. 


It  is  concluded  that  the  passive  acceptance  of  the  MOD  X TRUNCATION 
and  QUASI  STATIC  formulations  by  flight  control  analysts  should  be  justi- 
fied numerically  at  the  critical  flight  control  design  points.  In  addi- 
tion, it  is  concluded  that  the  full  implementation  of  the  integrated 
flight  control  benefits,  promised  by  the  CCV  concept,  has  a current  limi- 
tation due  to  the  possible  inappropriateness  of  the  "invacuum  normal  modes" 
as  elastic  coordinates  in  a highly  damped  system.  This  limitation  may  be 
removed  by  developing  a computational  method  for  -he  accurate  solution 
of  the  precise  EXACT  formulation. 
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Foreword 


This  Technical  Memorandum  was  prepared  under  work  unit  82191304.  The 
work  unit,  entitled,  "Analysis  Methods  for  Control  Configured  Military 
Vehicles",  is  divided  into  five  parts: 

1.  A literature  search. 

2.  A formulation  of  the  equations  of  motion  for  the  controls-free 
aircraft. 

3.  A simplification  of  the  equations  in  item  (2)  to  permit  a more 
rapid,  but  less  accurate  solution  for  flight  control  application. 

4.  Solutions  of  the  equations  in  items  (2)  and  (3)  using  USAF  air- 
craft as  test  cases  to  develop  criteria  for  the  selection  of  appropriate 
formulation  of  the  equations  of  motion. 

5.  A final  summary  report. 

The  work  reported  in  this  TM  summarizes  the  results  for  items  (1) , (2) , 
and  (3)  above. 

The  responsibility  for  the  accuracy  and  conclusions  presented  in 
this  Technical  Memorandum  rests  with  the  organization  that  prepared  it. 
This  Technical  Memorandum  has  been  reviewed  and  is  approved. 
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LIST  OF  MATRIX  SYMBOLS 


Matrix 

Symbol  Units 

A,  lbs/nrp 

AX  iWlTp 

Aj 


Definition 


Aerodynamic  forces  on  each  mass  due  to 
rigid  body  velocity,  "t/p  . 


Aerodynamic  forces  on  each  mass  due  to 
rigid  body  acceleration,  . 


Aerodynamic  forces  on  each  mass  due  to 
elastic  deformation  position. 


Aerodynamic  forces  on  each  mass  due  to 
elastic  deformation  rate. 


As  lbs/ Ip 


Aerodynamic  forces  on  each  mass  due  to 
elastic  deformation  acceleration. 


lbs/s; 


Aerodynamic  forces  on  each  mass  due  to 
position  of  ith  control  surface,  . 


lbs/£. 


Aerodynamic  forces  on  each  mass  due  to 
velocity  of  ith  control  surface, 


lbs/S; 


a Vlb 


Aerodynamic  forces  on  each  mass  due  to 
acceleration  of  ith  control  surface, 

• 

Flexibility  matrix  of  structure  canti- 
levered at  the  c.g.  of  the  aircraft. 


c 


Flexibility  matrix  of  the  structure  for 
a "free-free"  aircraft. 


ft 


D lbs/ip 

P lbs 

fc  lbs 


Displacement  of  mass  relative  to  its 
initial  position  as  measured  relative 
to  the  body-fixed  mean  axis. 

Structural  damping. 

Aerodynamic  forces  on  each  mass  due  to 
rigid  body  and  elastic  motion. 

Aerodynamic  forces  on  each  mass  due  to 
the  controls  &;(*)• 


Matrix 

Symbol 


Definition 


Body  axis  orientation. 

Residual  flexibility  static  aeroelastic 
correction  factor. 


ft/sec 


Generalized  coordinate  defined  as  in- 
vacuum normalized  mode  coordinate. 

Generalized  coordinate  defined  as  in- 
vacuum normalized  mode  coordinate  for 
r retained  modes. 

Generalized  coordinate  defined  as  in- 
vacuum normalized  mode  for  deleted 
modes. 

Body  axis  translation  rate. 

Combination  of  Y^and  "tr*  , ’"^pT- l^r'Top  J 


Rigid  body  mode  shape. 

Rigid  body  translation  mode  shape. 

Rigid  body  rotation  mode  shape. 

Mode  shapes  of  invacuum  normalized 
modes,  & - 4>,(t> x. 

Mode  shapes  of  retained  invacuum  normal- 
ized modes  r in  number. 


degrees 


Mode  shapes  of  deleted  invacuum  normal- 
ized modes. 

Motion  of  ith  control  surface. 
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1.0  INTRODUCTION 


Recently  the  flight  control  staffs  of  many  aerospace  research  and 
development  organizations  have  proposed  the  application  of  modern  flight 
control  concepts  early  in  the  design  cycle  of  elastic  aircraft.  The 
XB-70(-)  has  a gust  load  alleviation  system.  The  SST(2)  had  a "hardened" 
stability  augmentation  system  to  compensate  the  longitudinal  static 
stability  at  an  unstable  flight  condition.  The  B-l(3)  is  being  designed 
with  a "ride  quality"  augmentation  system  to  ensure  that  the  pilot  can 
effectively  monitor  and  control  the  B-l  when  flying  through  atmospheric 
turbulence.  Looking  to  the  future,  the  USAF  has  proposed (*0  the  Control 
Configured  Vehicle  (CCV)  concept  to  achieve  full  benefit  of  integrate  : 
flight  control  system  technology.  This  concept  requires,  among  other-  , 
that  (1)  the  reliability  of  a control  system  be  improved  to  a level  equal 
to  or  higher  than  the  reliability  of  the  primary  structure  and  (2;  the 
usually  adverse  effects  of  the  non-linearities  in  the  control  system  design 
data  be  overcome  using  high  authority  actuators  or  other  control  system 
hardware  and  software.  Then  the  problem  of  high  maneuver  loads,  low 
flutter  speeds,  static  instability,  etc.  on  a flight  vehicle,  that  is 
designed  to  be  aerodynamically  optimum,  can  be  minimized  or  removed  through 
the  use  of  a sophisticated  CCV-type  control  system.  Sone  examples  of  the 
CCV-type  control  systems  may  be  found  in  references  5,  6,  and  7. 


The  achievements  of  benefits  due  to  modern  flight  control  concepts 
depends  upon  an  accurate  mathematical  model  of  the  aircraft.  The  mathema- 
tical model  becomes  even  more  important  for  CCV  concepts.  This  is  because 
the  design  of  the  CCV  integrated  control  systems  requires  an  accurate  pre- 
diction of  the  rigid  body  and  structural  dynamics  of  the  aircraft  being 
stabilized  along  wxth  an  estimate  of  the  possible  errors  that  may  be 
included  in  the  dynamics.  In  general,  errors  in  the  dynamics  arise  from  two 
sources: 


1.  The  inappropriate  formulation  of  the  equations  of  motion. 

2.  The  ncn-precise  aerodynamic,  structural,  and  inertial  design  data 


set. 


It  has  been  traditional  to  attribute  most  analytical  errors  in  the 
flight  control  design  of  comtemporary  elastic  aircraft  to  the  design  data 
set  and  many  discussions  of  these  errors  are  contained  in  the  literature. 
Thus,  very  few  discussions  of  the  reasons  for  the  errors  in  equation  formu- 
lation are  to  be  found,  even  though  the  inappropriate  formulation  can  negate 
the  most  accurate  design  data  supplied  by  the  engineering  disciplines. 
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There  are  many  formulations  of  the  equations  of  motion  of  an  elastic 
aircraft  currently  used  in  the  aerospace  industry*  To  date,  most  of  the 
formulations  have  resulted  from  a priori  assumptions  or  have  been  forced 
upon  the  flight  control  design  due  to  expsdiencv.  These  types  of  inappro- 
priate formulations  result  in  moderate  to  large-sized  risk  factors 
in  the  flight  control  design.  The  risk  factors  are  manifested  as  the 
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large  gain  and  phase  margins  used  as  the  flight  control  system  design 
criteria  to  permit  large  adjustments  in  the  control  system  gains  during 
flight  tests  of  the  prototype.  More  complex  integrated  control  systems 
for  a highly  elastic,  high  performance  aircraft,  such  as  one  using  CCV 
concepts,  will  not  have  the  luxury  of  high  gain  and  phase  margins  be- 
cause of  their  complexity  and  their  use  of  each  aerodynamic  control  sur- 
face for  many  purposes. 

The  effect  of  an  inappropriate  equation  formulation  on  the  analysis 
of  the  integrated  flight  control  system  can  be  visualized  using  an  ex- 
ample. Figure  1.1,  adapted  from  reference  8,  portrays  a simple  feedback 
control  system:  the  schematic  is  presented  in  Figure  1.1a  and  the  sta- 

bility matrix  representation  of  the  schematic  in  Figure  1.1b.  In  both 
cases,  the  equations  of  motion  for  the  elastic  airplane  are  shaded.  As 
shown  in  the  schematic,  the  equations  of  motion  hold  a dominant  position 
in  the  design  of  integrated  control  systems.  Their  transfer  function 
relates  input  quantities  such  as  atmospheric  turbulence  and  gusts  and 
distributed  control  forces  and  moments  to  the  output  motion  variables 
that  determine  the  motion  of  the  elastic  aircraft  through  the  atmosphere. 
The  determination  of  the  numerical  value  of  all  the  elements  in  the 
schematic  may  be  achieved  by  the  evaluation  of  the  characteristic  equation 
of  the  determinant  of  the  stability  matrix  in  Figure  1.1b.  As  indicated, 
the  equations  of  motion  in  the  upper  lefthand  corner  of  the  matrix  in- 
fluence the  numerical  value  of  all  the  control  system  elements,  e.g.,  the 
filters,  forward  loop  compensation,  feedback  loop  compensation,  sensors, 
sensor  position,  etc. 

It  is  the  intent  of  this  paper  to  outline  the  source  of  those  er- 
rors in  the  integrated  control  system  due  to  the  formulation  of  the  equa- 
tions of  motion.  This  objective  is  accomplished  by  deriving  a unique 
and  precise  set  of  linear,  ordinary  differential  equations  of  motion 
for  a controlled  aircraft  flying  in  a gusting  or  turbulent  atmosphere 
(EXACT  formulation) . These  precise  equations  are  then  simplified  to 
five  other  formulations  using  identified  assumptions.  All  formulations 
reported  in  the  literature  and  used  within  the  aerospace  industry  can  be 
related  to  the  EXACT  formulation  or  to  the  five  simplified  formulations. 

A contrast  of  the  mathematical  model  of  each  formulation  to  the  actual 
physics  of  the  elastic  airplane  determines: 

• Which  formulation  is  the  most  appropriate  to  each  elastic  air- 
plane analysis. 

• What  qualitative  errors  are  included  in  the  integrated  control 
system  design,  if  an  inappropriate  formulation  is  selected  due  to 
cost  or  time  considerations. 

In  addition,  an  examination  of  each  of  the  equation  sets  indicates: 


O ^cv*eN\fcTic_  A Feets^>Acvc  Coktvcol  ^‘isthem 

Wtt  ={4 


/**  V • ' , '•,  \ 

/ /'L.v*  ' \ • 

/'  fc^OATKaMS  ' 

u 

-1/ 

r 

» 1 

y 

-w 

' op-  Motion 

P 

a 

OF  CorsTTKOLS,  FrSF  C 

* 

r 

d, 

Air chwpt 

/’  /'  .•* 

• 

' . ./  ] 
.«■’  . -•'  ’ * 

d- 

1 

(Qf. + l)£)  — ► •% 

Is  ! 

Sc 

>=< 

- — x * 

-1 

i«| 

FL  -l 

D. 

r 

* 

M -i 

*c 

c 

> 

K-\ 

1 

1 \ 

&<£ 

Oc.ommm<D 

K -t 

•eSp 

f 

O 

— 

A -J__ 

QaJ 

. * 

J 

b.  StAbM-iTy  M/wto*  op  a Svmplc  FeeoBACK  Con-ts-ov.  SvsreM 

^ 1,1  sSSfuF  E LAST  I (JT  A I K C RAFT  n^1' 1 ° ^ IN  ANALYSIS  0F  INTEGRATED  CONTROL 


t 


rr 


:‘-i 

;$ 


• What  problems  are  involved  in  implementing  the  formulations  on 
limited  size  flight  simulators. 

• What  formulation  is  particularly  appropriate  to  the  inclusion  of 
nonlinear  aerodynamics . 

1.1  History  and  Terminology 

The  analysis  of  the  dynamics  of  elastic  aircraft  is  traceable  to  the 
flutter  and  control  surface  divergence  aeroelastic.  problems  encountered 
on  the  ear^y  military  aircraft.  As  elastic  aircraft  became  more  complex, 
the  analysis  of  these  and  other  aeroelastic  problems  specialized  to  those 
engineering  disciplines  characterized  as  "flutter",  "static  structural 
loads",  "dynamic  structural  loads",  "configuration  aerodynamic  develop- 
ment", and  "aeroelastic  stability  and  control".  With  the  advent  of  more 
modern  aircraft  such  as  the  B-47,  B-52,  XB-70,  YF-12,  C5A,  747,  SST  and 
B-l,  the  interactions  of  these  specialized  disciplines  of  aeroelastic 
analysis  became  particularly  pronounced  during  the  design  of  the  flight 
control  system.  The  different  terminology  and  equation  formulations  of 
each  discipline  result  in  engineering  confusion  and  mis-natched  structural, 
aerodynamic,  and  stability  and  control  data. 

It  is  the  responsibility  of  the  flight  control  engineer  to  integrate 
the  data  created  by  the  aeroelastic  disciplines  to  create  a safe  and 
useful  aircraft  by  means  of  augmentation  systems.  These  aircraft  augmen- 
tation systems  are  used  to  improve: 

• Handling  qualities. 

• Ride  quality. 

• Static  stability. 

•Fatigue  life. 

• Flutter  margin. 

•Maneuver  loadability. 

•Atmospheric  gust  loadability. 

The  integration  of  the  diverse  aeroelastic  data  is  not  usually  performed 
by  any  of  the  other  disciplines  because  of  their  important  responsibilities 
to  the  detailed  structural,  aerodynamic,  or  stability  and  control  design 
of  the  aircraft.  An  additional  difficulty  encountered  by  the  other  dis- 
ciplines is  that  each  of  them  has  only  a partial  understanding  of  the  term- 
inology of  each  of  the  other  disciplines,  thus,  the  intercommunication 
required  for  integrated  control  system  development  is  difficult^). 
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Some  preliminary  formulations  of  the  equations  of  motion  of  an 
elastic  aircraft  have  been  completed.  Etkin(lO)  developed  the  equations 
of  motion  of  rigid  airplanes  and  outlined  a method  to  be  used  to  analyze 
elastic  aircraft.  More  specific  descriptions  of  Etkin’s  method  and  other 
methods  are  presented  in  references  11,  12,  13,  and  14  along  with  illus- 
trative applications  to  contemporary  aircraft.  There  are  some  publications 
discussing  the  equations  that  describe  the  motions  of  complex  elastic  air- 
craft configurations.  Unfortunately,  these  publications  specialize  the 
equations  to  an  undamped,  uncontrolled  elastic  aircra£t(l5)  or  to  a par- 
ticular formulation(16,  17,  18). 

Those  equation  formulations  commonly  found  in  the  literature  may  be 
categorized  as : 

QUASI  STATIC  - The  motions  of  the  structure  are  assumed  to  be  in  phase 
with  the  rigid  body  motions:  elastic  motion  acceleration  is 

instantaneous.  The  method  is  used  primarily  for  handling  qual- 
ity and  reduced  static  stability  control  system  design  for  elas- 
tic aircraft  with  wide  frequency  separation  between  the  rigid 
body  and  elastic  motions. 

EXACT  - The  motion  of  the  structure  is  determined  by  the  eigenvalue 
(root)  and  eigenvector  (mode  shape)  solutions  of  the  equations 
of  mot:' on  for  the  elastic  aircraft.  The  mode  shape  coordinates 
contain  complex  numbers-  The  accuracy  of  the  solution  is  limited 
by  the  existing  computerized  routines  that  calculate  the  complex 
number  eigenvalues  and  eigenvectors. 


MODAL  SUBSTITUTION  - The  motions  of  the  structure  are  assumed  to  be 
related  to  the  orthogonal,  invacuum  eigenvectors  (mode  shapes). 
All  eigenvectors  contain  only  real  numbers. 

RESIDUAL  STIFFNESS  - The  mode  shapes  representing  the  elastic  motion 
in  the  MODAL  SUBSTITUTION  formulation  are  separated  into  "re- 
tained" and  "deleted",  modes.  The  deleted  modes  are  represented 
in  the  dynamic  stability  analysis  as  quasi  static  aeroelastic 
corrections,  using  a correction  factor  related  to  the  deleted 
modes  and  the  stiffness  of  the  "free-free"  structure. 


RESIDUAL  FLEXIBILITY  - Similar  to  the  RESIDUAL  STIFFNESS  formulation, 
except  the  quasi  static  aeroelastic  correction  factor  is  related 
to  the  retained  modes  and  the  flexibility  of  the  free-free  struc- 
ture. 

MODAL  TRUNCATION  - The  deleted  modes  of  the  RESIDUAL  FLEXIBILITY  formu- 
lation are  not  represented  by  any  correction  factor.  This  is 
the  most  common  dynamic  aeroelastic  formulation  reported  in  the 
literature. 
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A unique  qualitative/semiquantitative  study  of  the  QUASI  STATIC 
RESIDUAL  FLEXIBILITY,  and  MODAL  TRUNCATION  formulations  was  conducted  by 
Macheai , Schwendler,  and  Pearce  for  the  USAF  and  reported  in  references 
16,  18,  and  19.  These  formulations  were  compared  to  the  RIGID  AIRPLANE 
and  EXACT  formulations  using  an  elastic  missile,  a low  aspect  ratio  air- 
craft, and  a B-47  aircraft  as  study  configurations.  The  aerodynamic  and 
structural  theory  applied  in  this  analysis  was  extremely  "crude",  but 
since  the  aerodynamic  model  was  held  constant  for  each  configuration,  its 
errors  are  hopefully  minimized  in  the  comparisons.  Typical  results  of 
the  study(^)  for  the  elastic  missile  anu  the  B-A7  aircraft  are  presented 
in  Figure  1.2,  a Bode  plot  of  the  log  amplitude  of  §/$c  versus  frequency. 
The  following  conclusions  were  noted  in  the  reference  16: 

• At  low  frequencies  (less  than  0.50  cps  in  the  Figure  1.2)  the 
QUASI  STATIC,  RESIDUAL  FLEXIBILITY,  and  MODAL  TRUNCATION  results 
are  nearly  coincident. 

*At  large  frequencies  (greater  than  0.50  cps  in  the  Figure  1.2) 
the  QUASI  STATIC  formulation  is  highly  inaccurate. 

•Neither  the  RESIDUAL  FLEXIBILITY  nor  the  MODAL  TRUNCATION  formu- 
lations accurately  approximate  the  EXACT  formulation  at  all  fre- 
quencies (above  0.50  cpc  in  the  Figure  1.2).  However,  the 
RESIDUAL  FLEXIBILITY  formulation  is  the  more  accurate  approximation. 

The  work  presented  in  the  reference  16  culminated  the  analytical  efforts  of 
the  USAF  to  represent  the  motion  of  elastic  aircraft.  Most  USAF  analytical 
work  in  this  field  was  discontinued  in  1962  due  to  the  absence  of  a digital 
computer  system  iarge  enough  to  solve  the  complex  equations  developed  using 
the  structural  and  aerodynamic  mathematical  models  that  describe  the 
physics  of  the  elastic  aircraft.  Instead,  the  USAF  concentrated  on  experi- 
mental programs  such  as  the  XB-70  GASDSAS  and  B-52  LAMS  to  demonstrate 
related  concepts. 

In  1965  - 1967  the  analytical  work  of  MacNeal,  et.  al.  was  studied  by 
NASA  and  the  Boeing  Company  and  reported  in  reference  12.  Here,  it  was 
recommended  that  the  RESIDUAL  FLEXIBILITY  formulation  be  implemented  using 
structural  data  from  the  structural  finite  element  program  such  as  NASTRAN, 
and  using  the  Woodward  aerodynamic  finite  element (20) . Also,  the  formu- 
lations of  reference  16  were  reworked  to  form  a amenable  to  the  large 
digital  computer  system  - the  CDC  6600  or  IBM  360.  The  resulting  program 
has  been  called  FLEXSTAB^l).  User  experience  with  the  QUASI  STATIC  and 
RESIDUAL  FLEXIBILITY  formulations  in  FLEXSTAB  has  not  yet  been  reported. 

1.2  Problem  Definition  and  Solution 

The  aerospace  industry  uses  many  of  the  formulations  described  in 
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CONFIGURATION  3,  HIGH  ASPECT  RATIO,  B-47  TYPE 


FREQUENCY  - Hz 


Section  1.1  to  design  and  develop  its  military  and  commercial  aircraft. 

Thus,  it  is  common  for  the  USAF  and  the  AFFDL  to  encounter  the  many  formu- 
lations (with  individual  company  variations)  during  the  proposal  evalua- 
tion and  the  contract  monitoring  phases  of  military  aircraft  development. 

As  examples : 

The  XB-70:  QUASI  STATIC  and  MODAL  TRUNCATION. 

The  B-52-LAMS  and  B-52-CCV:  MODAL  TRUNCATION. 

The  F-4  CCV,  Terrain  Follower,  and  Survivable  Flight  Controls: 

QUASI  STATIC  and  MODAL  TRUNCATION. 

The  YF-12  RIDE  QUALITY:  QUASI  STATIC  and  MODAL  TRUNCATION. 

The  F-lll:  QUASI  STATIC  and  MODAL  TRUNCATION. 

The  C-5A:  QUASI  STATIC  and  MODAL  TRUNCATION. 

The  B-l:  QUASI  STATIC  and  MODAL  TRUNCATION. 

The  F-15:  QUASI  STATIC  and  MODAL  TRUNCATION. 

\ 

The  proposed  SST  designs:  QUASI  STATIC,  RESIDUAL  STIFFNESS,  .j 

and  MODAL  TRUNCATION. 

There  has  been  no  previous  attempt  to  categorize  the  formulations  of  the  ' 

equations  of  motion  used  in  flight  control  system  designs  and  establish 
criteria  for  their  use.  As  a consequence,  the  limitations  of  each  formu- 
lation as  dictated  by  the  assumptions  required  to  derive  the  equations 
are  not  recognized. 

It  is  the  objective  of  this  Technical  Memo  to  classify  and  contrast 
these  formulations  relative  to  a unique  and  precise  formulation  of  the 
equations  developed  herein.  This  objective  is  accomplished  in  Sections 
2.0  through  5.0.  In  Section  2.1  and  2.2  the  nonlinear  and  linear  equations 
are  developed  to  describe  the  motion  of  a controls- fixed  elastic  aircraft 
relative  to  a non-inertial , body-fixed  coordinate  system.  The  linear  equa- 
tions are  modified  in  Section  2.3  to  include  the  effect  of  the  free  aero- 
dynamic control  surfaces  and  the  effects  of  atmospheric  gusts  and  turbulence 
on  the  elastic  aircraft;  the  resulting  equations  are  termed  the  EXACT 
formulation.  In  Section  3.0,  the  EXACT  formulation  is  modified  using 
listed  assumptions  to  simplify  the  mathematical  model  of  the  elastic  air- 
craft. These  simplifications  result  in  the  formulations  termed  QUASI 
STATIC,  MODAL  SUBSTITUTION,  RESIDUAL  STIFFNESS,  RESIDUAL  FLEXIBILITY,  and 
MODAL  TRUNCATION.  The  results  of  the  studv  are  presented  in  Section  4.0; 
the  conclusions  and  recommendations  of  the  study  are  summarized  in  Section 
5.0.  A List  of  References  is  contained  in  Section  6.0. 
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It  is  not  the  intent  of  this  paper  to  discuss  the  many  methods  avail- 
able within  the  government  and  industry  for  the  calculation  of  aerodynamic 
and  structural  influence  coefficients.  Thus,  it  is  assumed  throughout 
this  paper  that  methods  exist  to  calculate: 


• Distributed  aerodynamic  forces  on  the  surface  of  the  aircraft. 


• Mass  distribution  for  the  aircraft. 


• Stiffness  and  flexibility  matrices  for  the  free-free  and  canti- 
levered structure. 

•Other  data  as  required  for  the  flight  control  analyses. 

Also  the  effect  of  thrust  on  the  initial  shape  is  assumed  negligible. 

The  forces  due  to  thrust  perturbations  are  included  implicitly  in  the 
aerodynamic  force  terms. 

2.0  DERIVATION  OF  THE  EQUATIONS  OF  MOTION  OF  AN  ELASTIC  AIRCRAFT 

Before  the  flight  control  engineer  can  begin  development  of  the 
integrated  flight  control  system  of  an  elastic  airplane,  there  must  be 
an  accurate  description  of  the  aircraft  being  controlled.  This  descrip- 
tion is  termed  the  "equations  of  motion"  and  the  description  contains 
implicit  and  explicit  references  to  the  structural,  aerodynamic,  and  geo- 
metric properties  of  the  aircraft.  All  the  elastic  aircraft  equation  of 
motion  formulations  have  these  three  groups  of  data  in  common;  however, 
the  assembly  of  information  depends  upon  the  description  of  the  dynamics 
of  the  aeroelastic  problem. 

2.1  Selection  of  Motion  Coordinates 

There  are  many  descriptions  of  the  dynamics  of  elastic  aircraft 
discussed  in  the  literature.  These  descriptions  can  be  grouped  into  two 
categories  depending  upon  the  coordinates  used  to  describe  the  motions 
of  the  elastic  aircraft.  The  two  categories  of  displacement  coordinates 
are: 

•Displacement  coordinates  relative  to  an  inertial  axis,  i.e.,  an 
earth  fixed  axis  or  any  Galilean  related  axis. 

•Displacement  coordinates  relative  to  a non-inertial  axis,  i.e.,  an 
axis  "fixed"  to  the  accelerating  elastic  aircraft. 

These  two  axis  systems  are  shown  in  Figure  2.1. 
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2.1.1  Displacement  relative  to  inertial  coordinates 

In  the  first  analysis  grouping,  the  displacement  of  the  surface  of 
the  elastic  aircraft  is  measured  by  a vector  expressed  relative  to  an 
inertial  base  vector  system.  Typical  examples  of  this  approach  are  found 
in  references  Hand  17.  A flight  control  engineer  using  this  analysis 
approach  to  the  dynamics  requires  that  the  structural,  aerodynamic,  and 
geometric  vector  parameters  of  the  airplane  oe  expressed  in  the  inertial 
axis  base  vectors.  It  is  usual  in  most  of  the  literature  to  develop  the 
elastic  body  motions  as  a trajectory  of  the  body  and  the  surface  of  the 
body  relative  to  an  inertial  axis  translating  in  a straight  line  at  a 
uniform  velocity  relative  to  the  flat  earth.  The  uniform  velocity  and 
direction  of  motion  is  that  of  the  elastic  body  prior  to  the  onset  of  the 
disturbance  producing  the  elastic  motion.  This  analysis  method  is  typically 
used  by  those  specialities  concerned  with  '‘flutter"  and  "dynamic  struc- 
tural loads". 

2.1.2  Displacement  relative  to  non- inertial  coordinates 

In  the  second  analysis  grouping  the  displacement  of  the  surface  of 
the  elastic  aircraft  is  measured  by  a vector  expressed  relative  to  a 
non- inertial  axis  base  vector  system.  This  axis  system  experiences  both 
linear  and  angular  accelerations  relative  to  the  inertial  axis  discussed 
previously.  Typical  examples  of  this  approach  are  found  in  references 
12  and  15.  The  flight  control  engineer  utilizing  this  analysis  approach 
requires  that  the  structural,  aerodynamic,  and  geometric  vector  properties 
of  the  aircraft  be  expressed  in  the  body-fixed  base  vectors.  This  anal- 
ysis approach  is  typically  used  by  the  specialities  concerned  with  "aero- 
dynamic configuration  development"  and  "aeroelastic  stability  and  control". 

2.1.3  Contrast  of  inertial  and  non-inertial  coordinates 

The  flight  control  engineer  must  choose  the  better  of  these  two 
approaches  for  each  elastic  aircraft  flight  control  integration  problem. 

The  inertial  coordinate  analysis  is  a classical  mathematical  approach  and 
is  familiar  to  most  engineers  who  solve  structural  dynamics  problems.  On 
the  other  hand,  the  non-inertial  coordinate  analysis  possesses  consider- 
able flight  control  engineering  practicality: 

1.  Augmentation  system  criteria:  The  criteria  presently  used  to 

design  the  handling  quality  and  reduced  static  stability  augmentation  sys- 
tems are  expressed  in  terms  of  the  rigid  body  motions  and  aerodynamic 
stability  and  control  derivatives  measured  relative  to  the  non-inertial 
axis.  Since  criteria  for  the  integrated  flight  control  augmentation  sys- 
tem utilizing  CCV  concepts  are  not  yet  specified,  either  of  the  two  anal- 
ysis techniques  would  work.  However,  the  non-inertial  analysis  technique 
has  the  prerogative. 
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2.  Analysis  of  nonlinear  aerodynamic  effects:  The  integrated 

flight  control  system  design  points  usually  lie  at  the  extremes  of  the 
flight  placards  where  viscous  aerodynamic  effects  induce  pitch-up,  buffet, 
etc.  These  effects  are  currently  measured  in  the  non-inertial,  body-fixed 
axis  coordinates.  The  important  nonlinear  aerodynamic  and  control  system 
parameters  are  most  efficiently  programmed  in  terms  of  the  non-inertial 
coordinates,  thus  minimizing  the  required  computer  size,  computational 
frame  time,  and  required  axis  transformations. 


3.  Flight  simulator  analysis:  Pilot  work  load  and  disorientation 

are  very  important  considerations  for  the  design  of  an  integrated  flight 
control  system.  These  two  problems  are  usually  evaluated  in  fixed  base 
and  moving  base  flight  simulators.  Since  the  pilot  considers  himself 
a "mass"  attached  to  the  airframe,  he  evaluates  his  motions  relative  to 
some  body-fixed,  non-inertial  axis.  Usually  the  information  presented 
to  him  in  the  cockpit,  i.e.,  angle  of  attack,  sideslip  angle,  relative 
velocity,  et.,  are  directly  related  to  the  motion  of  the  non-inertial 
axis. 
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4.  Large  disturbance  maneuvers:  Large  disturbance  military  maneuvers  -S 
required  for  defensive  and  offensive  weapons  delivery,  are  very  often  de-  } 
sign  points  for  an  integrated  control  system.  These  large  disturbance  I 
motions  are  easily  described  in  terms  of  the  non-inertial  coordinates  and  | 
are  difficult  to  describe  in  terms  of  inertial  coordinates.  i 


Due  to  the  above  practical  considerations,  the  non-inertial  coordinate  I 
analysis  approach  is  adopted  in  this  Technical  Memo.  C viously,  both  j 
analysis  methods  must  give  the  same  results,  but  the  non-inertial  analy-  ; 
sis  is  the  analysis  method  most  easily  applied  to  the  flight  control  inte-  j 
gration  problem. 
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2.2  Controls-Fixed,  Elastic  Aircraft  Equations  of  Motion 

This  section  presents  the  equations  that  describe  the  motion  of 
an  elastic  aircraft  subject  to  aerodynamic  forces  that  are  independent 
of  control  surface  position  and  atmospheric  disturbances.  The  equations 
are  derived  for  the  case  of  the  aircraft  idealized  as  N "lumped  masses" 
each  having  3 translational  degrees  of  freedom  relative  to  a "body  axis 
system".  The  3N  rotational  degrees  of  freedom  may  easily  be  included  in 
this  formulation,  but  are  omitted  for  purposes  of  clarity.  The  body  axis 
system  is  defined  as  an  axis  system  "attached"  to  the  elastic  aircraft 
by  a means  to  be  specified  later  in  this  development.  Since  the  axis  is 
attached,  they  are  characterized  as  "non-inertial",  i.e.,  the  axis  experi- 
ences 3 translational  and  3 rotational  accelerations  whose  magnitude  are 
related  to: 

• The  elastic  motions  of  the  N masses  relative  to  the  body  axis. 

•The  aerodynamic  forces  on  the  system. 

•The  choice  of  the  body  axis. 

•The  location  and  orientation  of  the  body  axis  at  t = 0 (the  initial 
condition) . 

•The  gravitational  forces  acting  on  the  aircraft. 

•The  mass  distribution  of  the  elastic  aircraft. 


The  basis  of  this  derivation  is  Lagrange's  equation  rather  than  the 
Newtonian  equations  applied  by  references  12  and  15.  The  reason  for  this 
approach  is  its  simplicity  of  application  - a minimum  of  notation  and  of 
analysis  steps  are  required.  The  final  results  of  the  derivation  can  be 
checked  against  reference  12. 


The  Lagrangian  technique  to  be  followed  is  from  reference  22  and 
that  notation  is  used.  Recall  that  the  Lagrangian,  L,  must  be  written 
relative  to  an  inertial  axis  to  be  used  in  Lagrange's  equations: 


AGuA  _ ii. 


(2.1) 


where  is  the  Lagrangian  for  some  ith  mass  particle. 


(2.2) 
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^ is  an  appropriate  coordinate  vector, 
is  the  rate  of  change  of  L;  due  to  "q. 

-J^.(/v/)  is  the  time  rate  of  change  of  (^)  relative  to  an 
observer  in  inertial  space. 

^ is  the  impressed  force  vector  acting  on  the  system. 

"*S^j  is  the  angular  momentum  vector  of  the  system. 

» A 

^ is  a distance  from  the  axis  origin  to  a mass  point. 

Lagrange's  method  can  also  be  used  to  develop  the  equations  of 
motion  for  an  aircraft  whose  motions  are  measured  relative  to  a non- inertial 
axis,  i.e.,  base  vectors  J}B  =.  £&■/»(& and  • The  only 

additional  condition  to  be  applied  to  the  analysis  is  an  appropriate  def- 
inition of  (^): 


k^ 

time  rate  of 
change  relative 
to  inertial  axis 


time  rate  of 
change  relative 
to  non-inertial  axis 


to  X. 

time  rate  of 
change  of  base  vectors 


Consider  the  case  of  the  elastic  aircraft,  shown  in  Figure  2.1, 
which  is  experiencing  an  acceleration  through  a quiescent  atmosphere 
at  a translational  velocity  V«  and  a rotational  rate  J7.c Q • The  air- 
craft should  be  visualized  in  terms  of  a system  of  "lumped  masses"  inter- 
connected by  "springs"  and  excited  by  external  aerodynamic  and  gravita- 
tional forces. 


It  is  important  to  realize  that  the  body  axis  selection  for  the 
elastic  aircraft  case  is  very  difficult: 


• The  geometric  axis  attached  to  selected  masses  points  in  the  initial 
condition  is  non-orthogonal  at  subseouent  time,  i.e. , the  body 
reference  axis  system  (BS,  WBL,  WL)  becomes  non-orthogonal  for  the 
free  elastic  aircraft. 


•The  principle  axis  system  develops  translational  and  rotational  ac- 
celerations related  to  the  center  of  mass  motion  plus  the  elastic 
motion. 
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• The  attached  axis  system,  attached  to  a single  mass  such  as  an 
accelerometer  and  tangent  to  the  jig  body  axis  at  the  initial 
condition,  experiences  acceleration  proportional  to  the  center  of 
mass  motion  plus  the  motion  of  the  individual  mass  relative  to  the 
center  of  mass  motion. 

There  are  several  constraints  to  the  choice  of  the  body  axis  system: 

• The  aerodynamic  data  is  usually  available  only  in  the  jig  or  cruise 
geometric  axis  system,  the  stability  axis  system,  or  the  wind  axis 
system. 

•The  mass  distribution  and  geometry  are  specified  in  the  jig  or 
cruise  geometric  axis  system. 

•The  structural  properties,  stiffness  and  flexibility,  are  specified 
in  a structural  global  axis  system  that  may  be  coincident  with 
the  inertial  axis  or  the  non- inertial  axis  at  the  initial  conditions. 

It  is  assumed  for  purposes  of  this  derivation  of  the  equations  of  motion 
of  the  elastic  aircraft,  that  the  mass  distribution,  aircraft  geometry, 
and  structural  data  are  specified  in  the  cruise  body  axis  system.  Then, 
the  methods  of  reference  21  or  23  may  be  applied  to  transform  these  data 
to  a representation  of  the  aircraft  at  any  trimmed,  off-cruise,  M-q  flight 
point.  It  is  also  assumed  that  subsequent  structural  distortions  from 
the  trimmed  aircraft  flight  point  are  small,  i.e.,  changes  in  the  geometry 
of  the  planform  and  the  effects  of  material  plasticity  can  be  neglected. 

The  primary  structural  distortions  considered  are  those  contributing  sig- 
nificant forces  during  some  perturbation  motion  about  the  mean  motion  speci- 
fied by  the  initial  conditions  of  the  dynamics  problem. 

A more  microscopic  view  of  Figure  2.1  is  useful  in  defining  the  terms 
used  in  the  derivation  of  the  elastic  airplane  equations  of  motion.  One 
such  view  is  presented  in  Figure  2.2.  Here  the  inertial  axis,  the  body 
axis,  and  several  of  the  lumped  masses,  including  the  ith  mass  are  shown. 
F.ach  of  these  lumped  masses  represents  a portion  of  the  total  airplane  mass, 
i.e.,  fuel,  payload,  structure,  instruments,  etc.  Each  of  the  masses  has 
gravitational  and  aerodynamic  forces  acting  upon  it.  Each  of  the  masses 
has  accelerations  imposed  upon  it  proportional  to  the  net  system  accelera- 
tions (body  axis  accelerations). 

The  notation  used  in  Figure  2.2  is  defined: 

are  non-inertial  axis  attached  to  the  elastic  aircraft; 
base  vectors  and  . 

^ y £ are  the  inertial  axis  either  attached  to  the  flat  earth 
1 7 or  translating  rectilinearly  at  a uniform  velocity;  base 
vectors  X » ^ » and  >£  • 
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MICROSCOPIC  VIEW  OF  i th  MASS 
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m±s  a vector  distance  from  the  inertial  axis  origin  to  the 
body  axis  origin. 

p^(j0is  a vector  distance  from  the  body  axis  origin  to  the 
ith  mass. 

Y“.@t)is  the  vector  distance  from  the  inertial  axis  origin  to 
the  ith  mass. 

P^(t)is  the  vector  force  on  the  ith  mass. 

The  Lagrangian  for  any  mass  in  this  system  is: 


Li  = (%-n)  - 'A2  ^(px-fj) 


(2.3) 


where  VM;  is  the  mass  of  the  ith  "lumped  mass". 

K4“is  the  free-free  structural  stiffness  influence  coefficient  at 
the  ith  mass  due  to  a displacement  at  the  jth  mass  in  one  of  three 
possible  directions. 

Applying  the  definition  of  !*•  , shown  in  Figure  2.2, 


r;  - R-*p“  - Vca  4 4^ca  x Pc  (2,4) 


L;  ® Kij  (p;  - p^) 


(2.5) 


Equation  (2.5)  contains  6+3  N motion  variables. 


Substitution  of  L*  from  equation  (2.5)  into  equation  (2.1),  and 
letting ^ -a  , results  in  the  equations  of  motion  for  the  3 translational 
degrees  of  freedom  of  the  ith  mass  relative  to  the  non-inertial  axis: 


W;  (VCo>  2 J^ca*  T** 

♦ 

+iWx]p:  -v  Slc±»i5tco.xfl)) 


i K-iPi 


(2.6a) 
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For  axis  translation. 
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The  equations  describing  the  axis  rotation  are  more  difficult  to  develop 
since  the  axis  origin  is  not  attached  to  the  center  of  mass.  This  necessitates 
a careful  consideration  of  the  mathematical  operation  employed  in  equation 
(2.2).  Consider  first  the  angular  momentum  for  the  mass  system  in  Figures  2.1 
and  2.2  for  observers  stationed  in  both  the  inertial  axis  system  with  originOT, 
and  in  the  non- inertial  axis  system  with  originOp  : 

— N 

Kt  - 

Following  the  work  of  Goodman  and  Warner 
can  be  developed: 


o — w y*  • 

j k = $ f:  X m;b; 

T T 

the  relationship  of  rvrto  hp 


kf  = kt  - p * * M ? - ?xHr 


* 


where  p*  is  the  distance  from  oP  to  the  center  of  mass. 
M is  the  total  mass,  i.e. , M = ”2>  m;  . 


IT*  is  the  velocity  of  the  center  of  mass  relative  to  the  inertial 
axis. 


Vu  = M*-RXMa*-nf*x(£ 

< V «"  1 


where  Hj  are  the  moments  about  Oj  . 

are  the  moments  about  Op  . 

Ok*  is  the  acceleration  of  the  center  of  mass  relative  to  the 
inertial  axis. 


19 


-MfcM  taMiiii!atta^teaaMfflM^afeaiaAagMgaM!ii^BiMBaaBBaeiaMMaiMiMaa8iaaaliMMBttBaaaaMaiae 


i 


WSfF 

i 


^ r ?-*"  Wlf-1  ^ ra**vn!N*™rTO9fTr^  '---:™r'z^^^&^-'*-r 


% 


|j 


i 


e 


i 


if. 


i 


«« 

$> 


Thus,  ^f/ji-  7^  FTp  unless  the  non-inertial  axis  is  attached  to 
center  of  mass,pfc  = o>  or  R.  = o , or  p“*  is  parallel  to  R 


to  the 

, p"  = o.  or  k.  = o , or  p*  is  parallel  to  R . An 

alternate  form  for  ~np  is  that  presented  by  Milne^'  . Here, 


V2p*  - «u  i?i4i 


it* 


Kp  - h?  “ Max  * 2 


A convenient  redefinition  of  this  expression  for  Irtp  separates  the  agro- 
dynamic and  thrust  moments  flq,  , from  the  gravitational  moments,  , 

i.e. , 

4Ma 


h.=  Ma-  VcaXlm; 

1 4»l 


or,  alternatively, 
ha  = 


where  = Vc<v.  * Z 


" t ^/s,  + Vc„X  frn;  (PcaXft)- K 


for  an  axis  system  not  attached  to  the 
center  of  mass. 


ss  O for  the  axis  attached  to  the  center  of  mass. 

— % ■—  V 

In  terms  of  the  notation  in  this  paper,  hp  and  are  expressed  as: 

TL=Zft*  *>;&*+  + Jt<o.*p.) 

r ill 

M ^ ►»  H v 

= 2 f;  X "V  Vco.  +1  £ X Vs*  * i ft  * rv\;  (iW*^  ) 


1M 


Ha  - S^PA>*  4-jLaxfL  + h, 


N 


r%  f\  • H 

= iX(s^/st*V‘<0+  Z^3u^(p;xTta) 

i»«  -i‘s.1  4- 1 

1 421 

4 -I 

• 2 ^ m;  f . X (*pi4*  XJW)  + 


(2.6c) 


4*. 
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The  reader  should  observe,  at  this  point,  that  a simplification  of 
equation  (2.6c)  is^possible  by  noting  that  the  terms  f m,-  P‘/Si 
and  the  term  Wl;  will  cancel  one  another.  However,  this 

cancellation  will  Aiot  be  effected  to  measure  the  importance  of  the  attach- 
ment of  the  axis  orgin  to  the  center  of  mass  to  occur  later  in  this  paper 
in  section  2.2.2. 


2.2.1  Large  disturbance  equations  of  motion 

Equations  (2.6)  are  the  equations  describing  the  general  motion  of 
elastic  flight  vehicles  in  the  quiescent  atmosphere.  In  practice,  these 
equations  are  modified  and  simplified  to  facilitate  their  use  in  the 
analysis  and  synthesis  of  flight  control  systems.  The  description  of  the 
large  disturbance  motions,  to  be  defined  the  "large  disturbance  equations", 
is  required  for  two  purposes  in  these  analysis.  They  are  used  (1)  to  develop 
the  linear  equations  describing  the  small  disturbance  motions  and  (2)  to 
evaluate  in  detail  the  control  systems  designed  using  the  more  simple  and 
less  accurate  linear  equations. 

To  develop  the  large  disturbance  equations,  it  is  assumed  that  all 
the  time  dependent  motion  variables  can  be  separated  into  an  initial  value, 
plus  a perturbation  variable.  In  the  case  of  equations  (2.6a),  (2.6b),  and 
(2.6c),  this  separation  has  the  form: 


A substitution  of  the  perturbation  expressions  defined  in  equations 

(2.7) ,  into  the  motions,  related  by  equations  (2.6),  results  in  equations 

(2.8) .  In  equations  (2.8)  vector  products,  i.e.  ,o.xTd  , have  been  converted 
to  an  equivalent  indicial  -notation  form,  i>e-’ 

Axis  Translation: 

4 eJimn  (-flwEr*;  ni  + <**£*i  + 

where  i ranges  from  1 to  N;  where  j,  k,  1,  m,  n range 
from  1 to  3.  N 


(2.7a) 


(2.7b) 


(2.7c) 
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The  term  ea|  is  the  vector  permutation  symbol  commonly  employed  in 
the  indicial  notation  to  represent  vector  cross  products. 


Axis  Rotation: 


(2.8b) 


where 


(Ap.+Vip.) 


4 m;  e jj*  (3^+^)  e^mn  (P**; 

4 WV^j/4  ^Xi+PA;)  4 pi,;) 

- rr\;  <£^1  Sh,n  ( J>,^Om)(Piu;4 Pru) 

4 ftl;  p*t-  X'^4  4*p/4‘) 

- 3/nj  (?A,'+pA;^  (£*<'  +fww)^n'+i£0 


4 m;  eiM  e^n  (^^JCF^-vp,,;) 

-K-^) 


where  S*n  is  the  KronecVer  delta,  function,  where  i ranges  from  1 to  N; 
where  j,  k,  1,  m,  n ranges  over  1,  2,  and  3. 
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Elastic  Deformation  at  ith  Mass: 


+ X (%♦«*)($.+ + e-u 

N 3 

(2.8c) 

where  is  the  force  on  the  ith  mass  in  the  jth  direction  due 

to  a deformation  at  the  kth  mass  in  the  1th  direction,  and  where  no 
summation  on  i is  intended;  where  j,  k,  1,  m,  n range  over  1,  2,  and 
3 unless  otherwise  noted. 

In  equations  (2.8),  the  initial  values  of  the  parameters  are  F^*  } ; Ma  • Mg' 

M*-  ’ Fii  ’£);  « V-  , and  Jlj  ; the  perturbation  values  of  the  parameters  * 1 ^ 

are1^  M*.  , Mg ■ > M4- ,«f  , pjC  >~irj  and  **>  j . 

Equations  (2.8)  contain  both  initial  condition  motions  and  perturbation 
motions  from  the  initial  conditions.  It  is  of  advantage  in  later  analysis 
to  now  separate  the  two  motions.  This  can  be  accomplished  by  setting  all 
the  perturbation  parameters  to  zero  in  equations  (2.8): 

Axis  Translation: 


Axis  Rotation: 


where 


(2.9a) 


(2.9b) 


• * 

Up.  - wv  e^f  -vr^e^ 

+ ?i:  e^P*;  Jln 


£>i\t  MVCw!  &.1.V  W'AUU/.'i’.S  JttWiftMriS:  Vlrff  \<  '.tlW»K  .* 


u>.ikrt.Kt t±nr-*T>. 


2k 


J.......-~^»A'^tf<iri>Miai^ia^a.  tejjgl 


nr'-vr*"  '•"•KeKW'.TeR**: <7'™Z7*  • Jj» 


s;  r-r«?  WTren#" 


where 


hp.=  W;e.^t  (E^-t-j^V* 4tr<pA;)  (P&i^  +Pa;^4  p*;“£fc) 

4m*  <^k*  %efo,n  P«; Vn  '4  m»  ^jfc«  ^W*"  * P"^*  ffvO 

4W«*^ke  Pi«  Qinn^^  + W; £}  £{  (?A;4f&  )<%,nn (^"* P«;  +®*» £«* pd 
" CjjU  -ty* B;  Sbm  pn;  + ^mPni 

- wlx*  p/;  4-^A.Pf**  &mr\  Glm-Ru  +-P"'Pn«’  ^mP/i;4'(i)mfyi') 

4 m;eiA|  fei;  ft.  +piiPf;  +p^^y) 

“ p^ejemnEn;-^^ 

- Afh;e^^  (P/L  +pAi)  (Fm^n-V  p*.;-ftn  + pv»u'  J 

tV^p^-rOr^Eo-  *\rAp/;) 

4 rv>;  ^mPru' 

m>\ 

— • 


Elastic  Deformation  at  ith  Mass:  | 

■f- i ■=  y*;  -t-  pV  | 

(-typ/,  £/;  +<°i.  Pi;  ) ] 

* ^APA* 4 +^A  ffe)  + e^,n  (xpni  1 

4 rw  Pn;~] 

’< 

4 


Mr, 
^i=i 


jiipa 


(2.10c) 


The  large  disturbance  perturbation  equations  contain  reference  to  both  the 
initial  conditions  and  to  the  perturbation  motions.  Usually,  these  initial 
conditions  are  specified  or  determined  from  equations  (2.9)  in  each  problem 
of  interest  to  the  flight  control  engineer. 


2.2.2  Small  distrubance  equations  of  motion 

In  many  cases  of  practical  importance,  the  large  disturbance  equations 
are  simplified  to  a small  disturbance  form.  This  simplification  accom- 
polished  by  assuming  that  all  perturbation  variables  are  of  order  £ ,0(E.), 
where  £ is  a very  small  number.  Then,  products  of  the  perturbation 
variables  are  of  even  smaller  order,  i.e.,  O (£.  ) multiplied  by  O ( fc.  ) is 
£> ) , where  0(  £»  <<■  O ( £. ) . If  all  terms  of  0(  £,x)  or  smaller 

are  eliminated,  equations  (2.10)  are  linearized,  and  the  small  disturbance 
equations  result: 

Axis  Translation: 

f . = H Lvr-  co^Vt)]  * m-  fj;  ^kt 

+ fh;  (5*  -4  t0^  P(4*)  4 ^ A ( A.  fV>i  4 ^»wPrv,' ^ 
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Axis  Rotation: 


S = *ri  + % 


(2.11b) 


where 


Vy-vvi; e^{  (P^-*- p#;V«)  -»-  w;ejtt 

•VMi  Vn  -*-  i-  Pnj  Vy^ ^ 

+ to;ejkLpkeW^T*P^+,*-‘  e^^Aie^n 

- ^ejU^Eer s mn 

- ft);  Co p^'f  ) § n»n-^**  In;  +^-r  ej|?J!  jV**  + fi;^/) 


- ,3 rru p^.  e/p%tx p^;  J1  n 

- Xm;  ej  w e/ 

jvin.  (PrH;«V+  p«;-H  J 


= m;  e^$  ^L^.-tU^O  + m;  e^-i^e^JuE,; 


-Me 


Elastic  Deformation  at  ith  Mass; 


-Tjl  rr  rr\;  +^n  ^aVa)  -*-  fy  -U ej*e  ^fu  +(0a^ 

+ + ^kt$ke$nn (fl**p*;+<ArS*: ) 

^ ^ hi.  ^mrv  **  "Pru  ^ 

(2.11c) 


FVv?r?  7^*?  <7  =*  *.■  tZfrypz  mi*#! 


The  most  common  flight  control  problem  analyzed  using  the  small  disturbance 
equations  is  that  of  wings  level,  constant  velocity,  rectilinear  flight 
parailei  to  the  surface  of  the  "flat"  earth.  As  a consequence,  the  remainder 
of  this  paper  will  consider  that  analysis  problem.  For  these  initial 
conditions,  equations  (2.11)  are  greatly  simplified: 


Axis  Translation: 


^ M + ^^VjL  ^ ^ PV  + ejU  <^k  ZjLi  (2.12a) 


Axis  Rotation: 


M«.  - S>j  + M3i 


(2.12b) 


Kd  - ^ e<) M pA; V*  4 Vyi; <£j H +^e.w^ew?ft;Vt 

EV»*  + "U  *p\; 


Elastic  Deformation  at  ith  Mass: 


f.;  =,  m;  ai2c) 


Equations  (2.12)  are  expanded  in  equations  (2.13)  using  the  following 
definitions: 


"tr*  s J*-b 

2o  - 40,^+ 

s P **  Hf*  + 

fi  = P»; 


= + rX 


(a+S:)  - W+ZW;I^  = (Ai  **,)**  +(A«>0£fc 

“ te,+l>.)AB  -V  (Bl+U)^fc  4- 

(^*+c^  4-^; pi  = (C,-k:,)2$  ■vfe.+O^*  *4- 


where  u,  v,  w,  p,  q,  r;  Uo,  Vo,  Wo,  P,  0,  R are  the  standard  notation,  such 
as  defined  in  reference  10. 


Rigid  Body  Translation; 

ft -Mg  (©cos  ©J  - M(a+^W6-V'Vj  4 

fz4Mg$CAS©0  =r  M(v+tU6-pWft)  4 

4?3  - Mg  ( © 5>n  ©„)  = M (^-vpV0  -q\jb ) + 


a. 


a, 


+ -fCA 
4ve,-pC3 

*Pe--?c» 


(2.13a) 


Rigid  Body  Rotation 


Mai=r 

Wjj^+uS  ^j.-\^1s3*‘VCj 

• 

+ j>  ^ -Ixj  f 4 

2^*5;  *p3; 

^(v'.c,-u>cj-r(ut,e2-w„0 

u„^+uc3-w(i>,-^e,  ->m3j 

-Iycp+Ij^  -I^f  4 

»=|  1 <»  • 

♦ v(tf4Ct -V.C,)  -p(v.c,-tiC' 

f 

M.»  = 

! 

^ot>,  4lrC,-Udi)x-uCt  + 

- Vp-T3vi  +i»f-  + 

M M 

+p(u.?J-w.e,)-|(w.ct-v4) 

(2.13b) 

Elastic  Deformation  at  ith  Mass: 


?ti  - rn^te <«0,)  = h;(i^W.-rV,*jP3i.f £,.)  4^'p,;  +^K<llt#  Pll 


A*i 

•F^+mij^ss.  = M;  (irtni.-pVf  B„-  - pP3>.)  4 +t|K;lW  p,k 

H 


= ft;  (-Vi  4f V.  U.  4 P ) -U» ; p3;  4$RijUpiL 


(2.13o.) 


where  the  product  p^  is  the  stiffness  force  on  the  ith  mass  in 
the  l%  direction  (subscript  1)  due  to  the  vector  displacement, 
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of  the  kth  mass.  Similarly  for  and 

The  gravitational  forces  acting  on  the  system  appear  explicitly  in  equa- 
tions (2.13a)  and  (2.13c)  and  remain  implicitly  represented  in  equations 
(2.13c).  Following  the  convention  of  most  aeroelastic  analysts,  the  terms 
"axis  translation"  and  "axis  rotation"  are  replaced  by  "rigid  body  transla- 
tion" and  "rigid  body  rotation".  This  less  precise  terminology  has  become 
standard  even  though  the  term  "rigid"  is  certainly  misleading  in  elastic 
aircraft  analysis.  Its  use  illustrates  a commonly  encountered  problem  - 
often  one  aeroelastic  engineering  discipline,  e.g.,  flutter,  loads,  sta- 
bility and  control,  etc.,  generates  a unique  terminology  that  appears  in- 
consistent or  non-precise  to  other  engineering  disciplines. 

Up  to  this  point  in  the  development,  the  body  fixed  axis  system  has 
been  "fixed"  arbitrarily  to  the  body,  since  the  selection  of  an  axis  sys- 
tem often  depends  upon  the  problem  of  interest.  However,  in  elastic  air- 
craft analysis,  an  axis  that  facilitates  the  solution  of  equations  (2.13) 
must  be  selected.  Consider  the  axis  attached  to  the  center  of  mass  of 
the  flight  vehicle.  In  this  case  the  terms  C^.  C^,  and  all  terms  in 

the  first  box  of  equations  (2.13b)  are  zero.  If  in  addition,  the  axis 
that  is  attached  to  the  center  of  mass  is  also  the  "mean  axis"  (15), 
those  terms  in  the  boxes  in  equations  (2.13)  are  eliminated.  The  result 
is  an  immense  simplification  of  equations  (2.13):  the  only  remaining 

coupling  between  the  axis  translation  and  rotation  equations  and  the  elastic 
deformation  equations  occurs  in  the  aerodynamic  terms  ■f'-  » P-  • , and 

M.;  • * * 

This  paper  will  select  the  mean  axis  for  subsequent  analysis.  Since 

the  initial  orientation  of  the  mean  axis  is  a user  specified,  this  paper 

will  select  a mean  axis  orientation  such  that  ■©„  is  zero  and  such  that 

the  products  of  inertia  I and  I are  zero.  Equations  (2.13)  become 

xy  zy 


- M (•v'-v-'rU„-pW6) 

- Mf'vCr  -VpVo-^o') 


(2.14a) 


Ma,  - p - Xxj  £ 

Maj.  — ^ 

Mm  ~ “TX3P  +£1%* 


(2.14b) 


~ VH;  {u  +^W6  - *V0  +^?a;  - >r  tu ) +*;  f, ; ^K',  *4  fkk. 

- Nl;  (ir+Y'Uo-  fW&  +rP»;  “ p^;)  +l^«  Pjk 

- A;  (^  4pV0-  3 U0+  fP^  Pii)  *pjx  P<* 

* KB  I 

(2.14c) 


The  kinematic  relationships  between  the  inertial  axis  and  mean  axis 
provides  a relationship  between  the  Euler  angles  Cf>  , © , and  ■')£"  and 
the  rotation  rates  p,  q,  and  r: 


• * _».«»  • ;.o 

© = q , 0 = p + r yrnQi  , = r ^e£&o 


The  combination  of  the  kinematic  relationships  and  the  equations  of  the 
dynamics  can  be  written  in  a compact  matrix  notation  as  in  equations  (2.15): 


(2.15a) 


(2.15b) 


(2.15c) 
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^ ^ A,Vp  4-  A^^/p  4-  Aj + A-4  ap  + As  ip 


(2.15d) 


where  A and  contain  those  aerodynamic  force  derivatives  propor- 
tional to  the  rigid  body  translation  and  rotation. 

A^,  A^,  A,.  contain  those  aerodynamic  force  derivatives  propor- 
tional to  the  elastic  deformations. 

Combining  equations  (2.15a)  and  (2.15b)  and  summarizing  equation  (2.15c): 


[j  +H,Up4  Mir0p)  = $Tf 


K 1 


p 

where 


m - Tm  o i 

L° 


(2.16a) 


(2.16b) 


Equations  (2.16)  and  (2.15d)  are  the  equations  of  motion  for  a con- 
trols fixed  elastic  aircraft.  They  describe  the  linear,  small  perturbation 
motion  about  a reference  state  (initial  condition)  of  wings  level,  recti- 
linear flight.  A comparison  of  equations  (2.16)  to  equations  (6.118)  and 
(6.146)  of  reference  12  indicates  an  exact  correspondence.  Thus  the 
Newtonian  approach  of  reference  12  and  the  Lagrangian  approach  of  this 
Memo  have  resulted  in  the  same  equations  of  motion.  It  should  be  noted 
that  this  mathematical  representation  of  the  aerodynamic  forces  due  to 
rigid  body  and  elastic  motions,  i.e.,  A^  , A£  > A^  , A^  , and 

A is  unique  to  reference  12  and  this  paper.  The  alternative  method 
often  found  in  the  literature  is  a combination  of  A.  with  A^  and  of 
A^  and  A^  with  A^  . This  combination  permits  the  use  of4 the  many 

unsteady  aerodynamic  theoretical  methods  currently  found  in  the  literature. 
The  disadvantage  of  the  combination  is  its  implicit  dependence  upon  the 
frequency  (or  time)  variable.  This  implicit  dependence  requires  a costly 
iterative  frequency  solution  for  the  aerodynamic  forces  due  to  the  rigid 


body  motions  "Up  and  the  elastic  motions  JL 


ld  prior  to  or  simultaneously 
with  the  determination  of  the  eigenvalues  of  tne  equations  of  motion  of  the 
elastic  aircraft.  An  iterative  frequency  solution  is  not  required  for 
the  aerodynamic  forces  due  to  the  ~1/p  and  described  in  this  paper. 

However,  this  computational  savings  currently  limits  the  validity  of  the 
solutions  of  this  paper  and  reference  12  to  the  "low  frequency"  eigen- 
values. The  upper  frequency  boundary  (Figure  2.3)  of  the  validity  of 
the  "low  frequency"  approximation  has  not  been  determined,  but  from  the 
practical  viewpoint,  it  is  probably  sufficient  to  do  most  current  integrated 
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control  system  designs. 


Figure  2.3  Comparison  of  Low  Frequency  Approximation  to  Exact,  Unsteady 
Aerodynamic  Theory 


2.3  Control s-Free , Elastic  Aircraft  Equations  of  Motion  - EXACT 

Formulation 

A controls-free,  elastic  aircraft  analysis  problem  is  characterized 
by  active  aerodynamic  control  surface  motion  in  response  to  pilot  or  aug- 
mentation system  commands.  Two  common  reasons  for  these  necessary  com- 
mands are  atmospheric  gusts  and  atmospheric  turbulence.  The  equations 
for  the  controls-free  analysis  are  not  as  determinable  as  equations 
(2.16)  and  (2.15d)  for  the  controls-fixed  analysis.  The  reason  for  this 
is  the  dependence  of  the  turbulence  and  gust  aircraft  responses  to  the 
geometry  and  the  flight  Mach  number  of  the  aircraft.  This  difficulty 
can  be  overcome  by  assuming  a general  form  for  the  additional  terms  to 
be  added  to  equations  (2.15d).  The  net  results,  equations  (2.18),  or  (2.19), 
are  termed  the  EXACT  formulation. 

2.3.1  Structural  damping 

In  the  case  to  be  developed  in  this  section,  it  will  be  assumed  that 
structural  damping  is  viscous  in  nature  and  proportional  to  elastic  dis- 
placement rate,  . The  damping  force  is  represented  as  Dip  , 

where  D is  a damping  matrix  and,  in  general,  is  a dense  array.  It 
should  be  noted  that  the  "usual  assumption"  that  0 is  proportional  to 
frl,’  or  }<  "has  not  been  made,  thus,  the  effect  of  these  "usual  assump- 
tion" can  then  be  identified  in  the  development  of  the  special  formulations 
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in  Section  3.0.  In  the  event  that  0 is  proportional  to  fW4*  or  K 
the  following  identies  are  valid: 


m 

£?. 

i 


(faT  D (p  — o 


(2.17a) 


<#td  <t>  (2*17b) 

(2>T05  — O (2.17c) 


2.3.2  Atmospheric  gusts  and  turbulence 


The  representation  of  the  atmospheric  gust  or  the  atmospheric  turbu- 
lence often  depends  upon  the  configuration  being  analyzed.  Commercial  and 
military  aircraft  are  designed  using  different  gust  and  turbulence  mathe- 
matical models.  A given  military  aircraft  may  be  designed  to  satisfy 
several  gust  and  turbulence  mathematical  models,  depending  upon  the  extent 
of  its  flight  profile.  It  will  be  assumed  that  both  the  gust  and  turbu- 
lence can  be  represented  as  3 components  of  force  on  each  of  the  N masses. 
Also,  it  will  be  assumed  that  gust  forces  "fk  and  turbulence  forces 
have  the  same  element  ordering  as  aerodynamic  forces  p . The 
total  moments  and  forces  applied  to  the  rigid  body  are  then: 


(2.17d) 


The  general  form  of  and  allows  several  different  gust  forms, 

i.e.,  ramps,  steps,  1-^cos,  sin,  etc.  and  several  different  turbulence  forms, 
i.e.,  1,  2,  or  3 component  Von  Karman  or  Dryden  power  spectral  density 
representations.  It  is  only  required  that  the  spatially  defined  gust  and 
turbulence  be  transformed  to  the  body  fixed  mean  axis  prior  to  inclusion 
into  equations  (2.15d)  as  either  Pj  or  P*.  . Typical  transformation 

considerations  may  be  found  in  reference  (25). 
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2.3.3  Active  aerodynamic  control  surfaces  , | 1 

$ .1 

The  forces  due  to  specified  number  of  active  aerodynamic  control  j ' 

surfaces,  c in  number,  can  also  be  defined  in  a general  sense.  The  three  • 

components  of  control  surface  induced  forces  on  each  of  the  N masses  is  i i 

defined  to  be  . Again,  the  elements  of  are  arranged  in  the  ] 

same  order  as  those  in  f .It  will  be  assumed  that  the  elements  in  j 

•+*.  are  proportional  to  Sc  » S*.  , and  by  an  equation  to  be  ? 

developed  for  each  specialized  application,  i.e.,  ] 


KJj-  \ Ac; \ 4 \^t>\  +-A (2.17e) 


Inserting  the  general  representation  of  structural  damping,  atmos- 
pheric gust  and  turbulence,  and  active  aerodynamic  control  surfaces,  re- 
sults in  equations  (2.18): 


H ('Vp-i-  M.'ir p+  Mxf6pW  c t>r 


(2.18a) 


bip +-  K Jp  -tt;  (3p  4 0 Kp  -»■  H.'Vp 4 Mi_r0p  )) 

4 f i*  "f<-  4 

where  - A.Vp*  Af  4 A^lp  + As^f 


(2.18b) 

(2.18c) 


Alternatively: 
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where 


(2.19a) 


^4  - £ 

■ =A,-v,^AJ-Cf+/iJJ!f+A4if+A3*i!'p4-'Pc'vr3+Ft:  (2-l9b> 
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Equations  (2.19)  are  defined  to  be  an  EXACT  formulation  of  the  equations 
of  motion  of  a controls-free,  elastic  aircraft.  Inspection  of  these  equa- 
tions indicates  their  form  to  be  linear  ordinary  differential  equations 
requiring  a simultaneous  solutions  for  the  3N+6  unknowns  and 

subject  to  the  6 constraints  imposed  by  mean  axis  attached  to  the  center 
of  mass  of  the  aircraft. 

It  is  particularly  important  to  note  that  the  solution  of  equations 
(2.19)  will  result  in  both  complex  number  eigenvalues  and  eigenvectors: 

• The  real  and  imaginary  portions  of  the  eigenvalues  are  a measure 
of  the  stability,  circular  frequency,  and  phase  difference  of  the 
eigenvectors . 

•The  real  and  imaginary  portions  of  the  eigenvectors  reflect  the 
spatial  orientation  of  the  masses  associated  with  an  eigenvalue, 
i.e.,  the  masses  of  the  ith  mode  shapes  do  not  have  precisely  0° 
or  180°  phase  difference  in  position  when  the  ith  eigenvector  is 
complex  (reference  26  and  27) . 

It  will  be  shown  in  the  next  section  that  this  result  from  equations  (2.19) 
is  different  than  that  currently  calculated  by  the  "flutter"  and  "dynamic 
load"  engineering  disciplines.  If  the  flight  control  engineer  chooses  the 
typical  "flutter"  equations  of  motion  in  which  the  real  number  mode  shapes 
are  used,  the  phase  and  gain  relationships  of  any  control  system  that  is 
synthesized  may  be  in  error.  The  magnitude  of  this  error  depends  upon 
the  relative  magnitude  of  the  aerodynamic  terms  and  structural  damping 
terms  due  to  the  elastic  motion. 
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3.0  APPROXIMATE  FORMULATIONS 


The  solution  of  the  EXACT  formulation  requires  a large  digital  com- 
puter. A key  element  in  the  computerized  solution  is  the  routine  that 
determines  the  complex  number  eigenvalues  and  eigenvectors.  Presently, 
these  routines  can  calculate  a "small  number"  of  complex  number  eigen- 
values and  eigenvectors  with  a high  degree  of  accuracy.  ' As  the  number 
of  complex  number  eigenvalues  increases,  the  accuracy  of  the  routines 
deteriorates  due  to  the  limited  number  of  significant  digits  that  the 
computer  can  retain.  In  the  case  of  a solution  to  the  EXACT  formulation, 
even  those  calculations  using  double  precision  will  eventually  suffer, 
since  a typical  value  of  an  N is  300  or  approximately  900  eigenvalues. 
Consequently,  an  accurate  solution  of  the  EXACT  formulation  of  the  equa- 
tions of  motion  is  unrealizable  for  all  but  simplified  analysis  problems. 

In  this  section,  five  approximations  to  the  EXACT  formulation  of 
the  equations  of  motion  for  the  controls-free  aircraft  are  constructed 
to  permit  a cheaper  and  faster  determination  of  the  unknown  motions. 

These  approximations  are  developed  by  reducing  the  number  of  unknowns  and 
equations  through  the  use  of  assumptions  that  may  be  valid  only  for  a 
specific  elastic  airplane  problem.  In  a mathematical  sense,  there  is  no 
single  set  of  these  approximate  equations  that  are  valid  for  all  aircraft 
or  a single  aircraft  at  all  flight  conditions.  In  an  engineering  sense, 
one  or  more  approximate  sets  of  equations  may  be  valid.  In  each  of  the 
five  formulations,  a matrix  equation  similar  to  the  EXACT  formulation, 
equation  (2.19),  will  be  developed  as  the  final  result. 

The  reduction  in  the  number  of  the  equations  and  unknowns  often  re- 
sults in  equations  of  motion  that  are  from  10  to  50  in  number.  This  order 
of  reduction,  e.g.,  from  100's  to  10's  in  number  is  sufficient  to  over- 
come the  numerical  difficulty  in  the  digital  computer  routines  used  to 
calculate  the  complex  number  eigenvalues  and  eigenvectors.  If  the  reduc- 
tion in  equation  number  and  complexity  is  large  and  excess  central  memory 
storage  space  becomes  available,  the  nonlinear  aerodynamic  and  structural 
data  may  be  incorporated  into  the  analysis.  These  nonlinear  terms  are 
not  considered  explicitly  in  equations  (2.19),  but  may  be  easily  included 

There  are  two  means  to  effect  the  simplification  of  the  EXACT  formu- 
lation: 


1.  The  eigenvalue  solution  of  the  EXACT  formulation  equation  can  be 
truncated  at  some  selected  frequency.  The  effects  of  the  neglected  higher 
frequency  eigenvalues  must  then  be  "filtered"  from  the  integrated  control 
system. 


2.  The  EXACT  formulation  equations  can  be  simplified  to  eliminate 
the  complexity  by  reducing  the  number  of  equations  and  unknowns,  prior  to 
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their  solution. 


There  is  no  engineering  experience  at  the  present  tine  that  would  indicate 
which  of  the  two  simplifications  is  the  more  appropriate  for  elastic  ve- 
hicles at  all  M-q  flight  conditions. 

The  second  simplification  is  that  most  commonly  employed  and  reported 
in  the  literature.  When  this  simplification  is  selected,  the  flight  con- 
trol engineer  must  decide  upon  a static  or  dynamic  aeroelastic  formulation 
appropriate  to  his  elastic  aircraft.  This  decision  requires  a certain 
amount  of  engineering  insight  based  upon  experimentation  and  previous  ex- 
perience. It  is  the  second  simplification  that  will  be  discussed  for  the 
remainder  of  this  paper;  the  first  simplification  is  currently  being 
investigated  by  the  AFFDL/FGC. 

The  reduced  set  of  equations  to  be  considered  in  this  Section  are: 

• QUASI  STATIC 

• MODAL  SUBSTITUTION 

• RESIDUAL  STIFFNESS 

• RESIDUAL  FLEXIBILITY 

• MODAL  TRUNCATION 

3.1  QUASI  STATIC  Formulation 

The  QUASI  STATIC  formulation  is  the  classical  method  used  by  stability 
and  control  and  flight  control  engineers  prior  to  the  advent  of  highly 
flexible  aircraft.  This  formulation  is  the  most  rapid  and  the  cheapest 
to  solve.  Due  to  the  small  number  of  equations,  it  is  ideally  suited  to 
a nonlinear  aerodynamics  analysis  based  upon  experimentally  measured  data. 
An  illustration  of  this  approach  is  found  in  reference  12. 

The  derivation  of  the  QUASI  STATIC  equations  begin  with  equations 
(2.19): 
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(2.19b) 
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where 


^ s At'”'Jp'»  Ai/l5p  *f  A3  j)p  4 A4^p  4 As3lp  4-  ^c.4^  4<V 


These  equations  are  simplified  with  assumptions  A1  and  A2 

A1  the  aerodynamic  forces  proportional  to  cfip  and  <0^,  are  zero 
i.e.,  Mf  = Astp  =0.  F 

A2  the  structural  inertial  and  damping  forces  are  zero,  i.e., 

n»3TP  = Di)p  = 0. 

Equation  (2,19)  becomes: 
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Alternatively: 


(3.1a) 


4 w;2^rlp  + + rn;  5vrM,i  4*fc  Jp  ~ ft  (3.1b) 


where  now 


f}  - A,Trp4  t 


(3.1c) 


AO 


Substitute  ft  from  equation  (3.1c)  into  equation  (3.1b),  collect  the 
common  terms,  and  then  determine  <fip  : 

-*-M;f>M,/irp  +-  ^ K«fp  - ft  - 

“ A.TTpi- /liTTp  +/l3i)p 

«0p  - Tk-AaT'  [(Ai  - 4 (A,  - m^Mj'Vp 

“ + ^ (3.2) 


Substitute  of  equation  (3.2)  into  equation  (3.1c)  and  group  common 


terms: 


ft  - (A,  4 A^lK'Aj'KA.-m^M^lTp  + (Ac+A s(K-A3VW~mA-3>))lrp 
“ A^(K-A3)  + [z+A^IK-Ai) 

* o /I  «J\ 


(3.3) 


The  QUASI  STATIC  formulation  is  presented  in  equations  (3.4): 


ft  = U + AjfK'A^fA,-^  +A,/Vp  +£.4^4  ft) 

4 A3(V,-A3^  (-nv5-vf  -rv*;  $Mt.rftp 


(3.4a) 


(3.4b) 


E 


All  explicit  reference  to  elastic  motion  dp  has  been  eliminated 
in  equation  (3.4).  The  factor  (i+A*  (K-A3r')  is  a "quasi  static"  aero- 
elastic  correction  factor  applied  to  the  rigid  airplane  forces  and  moments. 
The  factors-fo^  » and  are  "inertial  relief"  corrections. 

An  alternate  formulation  exists  in  the  literature.  ' It  is  based  upon 
the  availability  of  a flexibility  matrix  for  the  cantilevered  structure 
and  the  construction  of  the  "free-free"  flexibility  matrix  6 . This 

formulation  is  repeated  below  from  reference  12.  Since  the  formulation 
in  reference  12  describes  motion  of  a controls-fixed  elastic  aircraft  in 
quiescent  air,  editorial  prerogative  has  been  assumed  and  the  terms 


have  been  added  to  the  equations. 
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(7  is  the  flexibility  matrix  for  the  cantilevered  structure. 

An  individual  applying  the  QUASI  STATIC  formulation  has  in  effect 
assumed  that  the  elastic  deformation  occurs  instantaneously  and  in  phase 
with  the  axis  system  motions.  This  representation  is  usually  valid  for 
the  design  of  handling  quality  and  reduced  static  stability  augmentation 
systems  of  elastic  aircraft  at  very  low  frequencies,  or  for  elastic  air- 
craft having  a large  frequency  separation  of  its  rigid  body  and  elastic 
deformation  motions.  The  formulation  is  not  generally,  valid  for  the 
design  of  the  other  CCV-type  control  systems. 

The  QUASI  STATIC  formulation  has  as  its  advantages: 

1.  The  solution  of  only  6 equations,  permitting  many  M-q  flight 
conditions  to  be  rapidly  analyzed. 

2.  The  inclusion  of  nonlinear  aerodynamic  effects  to  flight  controls 
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analyses  using  aeroelastic  correction  factors,  e.g. , '-*•«*  /Q*,'C',D, 
• «c- 


3.  The  convenient  implementation  on  limited  size  analog  or  digital 
flight  simulators. 


3.2  MODAL  SUBSTITUTION 


The  MODAL  SUBSTITUTION  formulation  is  classified  as  "dynamic  aero- 
elastic",  since  the  elastic  motion  is  not  assumed  to  be  in  phase  with  the 
rigid  body  motions.  It  is  assumed  that: 


A3  structural  damping  forces  are  negligible,  i.e. , very  small 
compared  to  the  structural  stiffness  and  structural  inertial 
forces. 


A4  aerodynamic  forces  due  to  the  elastic  deformation  are  negligible, 
i.e.,  small  compared  to  the  structural  stiffness  and  structural 
inertial  forces. 


Assumptions  A3  and  A4  permit  a coordinate  transformation  using 
the  invacuum,  orthogonal  modes  of  vibration.  This  substitution  results 
in  the  simplification  of  the  complex  number  eigenvector  solutions  of 
equations  (2.19)  to  real  number  eigenvectors.  The  procedure  for  accom- 
plishing this  transformation  is^ 


1.  Recall  equations  (2.18b): 


t)JL  + KcSp  — -W;  (^p  + Mi'Vp 


(2.18b) 


2.  Represent  dtp  as  an  ensemble  of  invacuum,  orthogonal  modes 
' , i.e.,  JL  = <$  ‘ — < - - - 


of  vibration,  i.e.,  tip  = U.  , where  = * 0. 

3.  Post  multiply  equations  (2.18b)  by  0r  . 


4.  Apply  assumptions  A3  and  A4 


The  result  is  presented  in  equations  (3.5): 

ortno4oMAuty 
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(3.6a) 


where 


<$Tm;  — 

<tT  K = i. 


(3.6b) 


5.  Determine  the  eigenvalues  and  eigenvectors  satisfying  equation 
(3.5)  such  that  (3.6a)  and  (3.6b)  are  valid,  i.e.,  find  the 
eigenvalues  of  the  determinate  3,'ry'  **• J 1 = 0*  This 

eigenvalue  problem  is  vrell  documented  in  reference  28. 


Of  the  , in  the  general  6D  case,  3N-6  in  number  are  real  non-zero 
numbers;  these  that  are  real  are  used  to  determine  the  value  of  0.  ^Six 
of  the  J\;  are  zero  due  to  the  use  of  the  stiffness  matrix  Ic  for  the  free- 
free"  structure.  Physically,  the  zero  eigenvalues  mean  that  the  elastic 
displacements  are  known  within  an  arbitrary  3 rotations  and  3 translations 
of  the  elastic  body  ensemble.  The  amplitudes  of  the  translations  and 
rotations  are  zero  by  definition  of  the  mean  axis.  The  degrees  of  freedom 
may  thus  be  removed  bv  combining  equation  (3.5)  with  equations  (2.18b)  and 
(2.18c): 


H 4 H.TTp 

(3.7a) 

Wi  ix  4-  A u. 

= Cf>T($ -c^D^u. 

(3.7b) 

~ j 

fttGp  + u + A4  $ + A<,4>tL 

(3.7c) 

Note  that  this  coordinate  transformation  of  Ap  to  has  reduced  the 
number  of  equations  and  unknowns. 

Equations  (3.7)  can  be  rewritten  into  the  matrix  format  similar  to 
equations  (2.19): 
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(3.8b) 
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The  equations  (3.8)  are  termed  the  MODAL  SUBSTITUTION  formulation. 

This  formulation  has  reduced  the  complexity  of  the  problem,  i.e.,  complex 
number  eigenvectors  in  equation  (2.19)  are  reduced  to  prescribed  real 
number  eigenvectors  in  (3.8).  However,  this  reduction  has  not  eliminated 
the  problem  of  determining  a large  number  of  complex  eigenvalues  and 
complex  number  eigenvector  magnitudes.  The  "inaccurate"  complex  number 
eigenvalue  routines  required  for  the  EXACT  formulation*  must  still  be  used. 


3.3  RESIDUAL  STIFFNESS  Formulation 


The  RESIDUAL  STIFFNESS  formulation  reduces  the  number  of  equations 
and  unknowns  that  are  associated  with  the  elastic  motion  in  the  MODAL 
SUBSTITUTION  formulation.  This  reduction  may  be  limited  to  only  one  mode 
shape  or  to  most  of  the  3N-6  mode  shapes.  The  reduction  is  accomplished 
by  the  engineer: 


1.  Noting  that  a large  separation  in  the  frequencies  of  the  mode 
shapes  has  occurred. 


2.  Deciding  that  the  problem  need  include  only  those  eigenvalues 
that  are  less  than  some  selected  frequency. 


The  RESIDUAL  STIFFNESS  formulation  then  represents  dynamically  all 
the  modes  retained;  it  "statically"  represents  all  the  modes  deleted. 


First,  separate  <*P  into  retained  and  deleted  modes: 

<$P  = td,  <&.]  + 

where  d,u,  are  modes  to  be  retained. 

are  modes  to  be  deleted. 

The  <$p  of  equations  (3.9)  is  substituted  into  equation  (3.8): 

‘M  o O o [MM,*,  o o Iff  7 G olO 

° “ p>  rU  0 !krH  6 ?M 

o o to,  o / u, 

.o  o o m,  \ut 


(3.9) 


EL  "* 


O U.7+  O w 0<*  /,*;  j o k o 

O mJVu, } L o c^TdqS  <&rD(*  1°  o *■' 


o o 


LO  o o ) L o <5^dqS  <arD*J  1 

where  W , •=  Cp? A,  - $ T JK 
= (tfmj  - C ^tk 


/0,TO 

l 

(3.10a) 


+ A3(<&U,+<&LlJ  + A+O&iUflWO 
4 As  ( Cb  u, + <A.ut)  * £c.  * 


(3.10b) 


In  order  to  reduce  equations  (3.10)  to  the  RESIDUAL  STIFFNESS 
formulation,  the  following  assumptions  are  required: 

A5  the  aerodynamic  forces  due  to  deleted  modal  velocity  and  accelera- 
tion are  zero,  i.e.,  i = A£C^yx  = 0. 

A6  the  structural  damping  and  inertial  forces  due  to  deleted  modal 
deformation  are  zero,  i.e.,  $,TD = 0,  = 0. 

A7  the  structural  damping  of  the  retained  modes  on  the  deleted  modes 
is  zero,  i.e. , $JD  d>,  = 0. 

Applying  the  assumptions  to  equations  (3.10)  reduces  them  to  the 
form  of  equations  (3.11): 
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4 ^c.  4 'Pq  + {qj- 

(3.11c) 

Equation  (3.11b)  is  solved  for  U.fc  and  the  result  is  substituted  into 
equations  (3.11c): 


-Tj  ^ u,  -v-/U<£,u.-+  As$,u, 


Equation  (3.12a)  has  a common  factor,  Tx  , appearing  on  both  the  right 
and  the  left  side  of  the  equality.  A grouping  of  this  factor  to  the  left 
side  of  the  equality  and  removing  the  multiplication  factor  by  inversion 
reduces  equation  (3.12a)  to  the  equation  (3.12b) 


a<kV<fe.T]  1 (^‘^p-tAaVp  + A4  <£,  u, 

4-  4 ^4^4^  ) 


(3.12b) 
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Summarizing  equation  (3.12b)  and  equation  (3.11a): 
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(3.13a) 
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tAs<AU,  ) (3.13b) 

Equations  (3.13)  are  the  RESIDUAL  STIFFNESS  formulation.  They  are 
6 + r in  number,  where  r is  the  number  of  retained  modes,  0,  . A 

substantial  simplification  of  the  MODAL  SUBSTITUTION  formulation  (equations 
(3.8))  has  been  accomplished.  The  3N  number  of  equations  and  unknowns 
has  been  reduced  to  6+r  number  of  equations  and  unknowns.  The  cost  of 
the  simplification  is  that  the  elastic  airplane  must  satisfy  assumptions 
A3  and  A4  as  in  the  MODAL  SUBSTITUTION,  and  in  addition,  must  satisfy 
assumptions  A5,  A6  , and  A7  . The  advantage  of  this  technique  is 
that  r can  be  selected  such  that  the  6+r  equations  are  less  than  50,  thus, 
the  "inaccuracy"  of  the  complex  number  eigenvalue  computer  routine  can 
be  minimized.  The  disadvantage  of  the  technique  is  that  all  the  mode 
shapes,  0,  and  <&>  must  be  calculated. 

3.4  RESIDUAL  FLEXIBILITY  Formulation 

The  RESIDUAL  FLEXIBILITY  formulation  eliminates  the  necessity  to 
calculate  0*  required  in  equation  (3.13)  of  the  RESIDUAL  STIFFNESS 
formulation.  The  elimination  is  accomplished  by  redeveloping  equations 
(2.18a)  using  the  "free-free"  flexibility  matrix.  An  excellent  descrip- 
tion of  this  formulation  is  found  in  references  12  and  16.  The  expres- 
sion equivalent  to  equation  (2.18a)  is  equation  (6.129)  of  reference  12: 

A?  - -<3  r*;  $*p-»  0 Ovpt  M.-Up*  Ma.r ' )j-  5 i'JP  + ah 

(6.129, 

r - , Ref.  12) 

where  C * l?  - <J>  M (P  H4;|C  is  the  "free-free"  flexibility  matrix  and 
Cf  is  the  flexibility  matrix  for  the  aircraft  structure  restrained  (at 
the  c.g.)  against  rigid  body  translation  and  rotation.  The  terms  CDtPp-K?-?; 
have  been  added  using  editorial  prerogative. 

The  removal  of  <k  k;'<£  from  equation  (3.12a)  is  accomplished  in 
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Collecting  terms  that  contain 

= Lc  - 4 h;1#]  (-*$  - & ^ f; ) 

- 5 r>.4  <£  ('VpHM.'Vp+Mxfbf  ) 


Defining  R.  = d - <$>, ){,  (f>, T : 

- C'm;^  (^M.^fp+M-tTop ) (3 . 16) 

4.  Substitute  Jp  = tf/Hi-t  (f>iM\.  from  equations  (3.9)  into  equa- 
tion (3.16)  and  recall  A5  , A6  , A7  2 


K^f+  Mx  Top] 


- R (- — Oq6(ut  +?; ) 

-C^$/^f**Mtnrp+^r4)  (3-17) 
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Substitute  equation  (3.18)  into  equation  (3.12a),  collect  the  terms 
on  left  side  of  equation,  and  remove  the  common  term  by  inversion: 

■f)  - A, -t  A\-v-p  -f  ^3^,  U,  -f^4du, 

+ Ai K(-r^'^u.  - D$,u,  +-T0 

- Abc  m;  ^ ('Uj.i  ) 


T“*  # | 

• x = j.1  - AjrJ  ( A/ u"p  4 Aj.'V'p  4 + ^4^^,  +As<$, u, 

^c.  (~*H«  ^ - P$. ) 

-Ai  ^ m„-  $ /'Vp  4 Mt-V-p*  H^Vapj)  (3.19a) 
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Equations  (3.19)  are  the  RESIDUAL  FLEXIBILITY  equations  that  describe 
the  motion  of  a controls-f ree  elastic  airplane.  As  in  the  RESIDUAL  STIFF- 
NESS formulation,  the  number  of  unknowns  is  6+r,  where  r is  the  number  of 
retained  modes  <j£,  . This  formulation  does  not  depend  upon  4.k?d>J 

as  did  the  RESIDUAL  STIFFNESS  formulation,  thus,  <f>t.  need  not  be  calcu- 
lated. The  penalty  for  using  the  RESIDUAL  FLEXIBILITY  formulations  are 
the  necessities  to: 

1.  Determine  C for  the  aircraft. 

2.  Calculate  more  numerous  matrix  products  in 

The  relative  cost  advantage  of  the  RESIDUAL  STIFFNESS  and  RESIDUAL 
FLEXIBILITY  formulations  is  unknown  at  this  time  and  is  being  investigated 
further. 

3.5  MODAL  TRUNCATION  Formulation 

The  most  common  formulation  used  to  solve  the  controls-f ree,  elastic 
aircraft,  dynamics  problem  is  MODAL  TRUNCATION.  This  formulation  is 
currently  applied  to  the  B-52  CCV  (reference  li)  and  the  B-l  (reference  29) 
aeroelastic  analyses.  This  formulation  reduces  the  computational  steps 
required  in  the  previous  formulations  by  assuming  A8  , and  A9  below, 
in  addition  to  A3  thru  A7  assumed  previously: 

A8  the  structural  spring  forces  due  to  U*.  are  zero,  i.e.,  Zt-t<4u=  0. 

A9  the  aerodynamic  forces  due  to  Uu  are  zero,  i.e.,  AjC^_ux=  0. 

The  application  of  assumptions  A3  through  A9  to  equations  (2.19) 
reduces  them  to  the  following  from: 
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A = A^p-f-  A^-lTp  -4  A s^U,  4 At^/U,  4 ASC^U,  + ft.-v£j+£*-  (3.20b)  j 

i 

Inspection  of  these  equations  indicates  all  reference  to  <^_<-U.  has 
been  deleted.  This  reduces  the  number  of  equations  and  unknowns  to 
6 + r , where  r is  the  number  of  retained  modes  cf>, 

This  formulation  has  the  advantage  of  a simplified  aerodynamic 
representation.  In  addition,  the  required  computer  size  and  the  computer 
central  processing  time  are  smaller  than  either  residual  formulation, 
thus  allowing  more  analysis  for  the  dollar. 

The  disadvantage  of  the  formulation  is  that  it  has  a narrower  range 
of  applicability  than  any  of  the  other  dynamic  aeroelastic  formulations. 

Thus,  its  use  in  the  design  of  integrated  flight  control  systems  should 
be  verified  with  quantitative  and  qualitative  analyses  that  test  the  as- 
sumptions A3  through  A9 
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4.0  RESULTS  I I 
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The  primary  results  of  this  study  are  the  six  formulations  of  the  i 

equations  that  describe  the  small  disturbance  motions.  These  six  formula-  1 

tions  are  summarized  in  Table  4.1.  Here  the  similarity  of  the  formula-  { 

tions  is  readily  apparent  as  well  as  the  reduction  in  the  number  of 
equations  that  occurs  with  each  of  the  successive  approximations.  As 
shown  in  the  table,  the  EXACT  formulation,  equations  (2.19)  are  3N+6  in 
number.  The  MODAL  SUBSTITUTION  formulation,  equations  (3.8),  reduces  the 
number  of  equations  to  3N  in  number,  while  the  RESIDUAL  STIFFNESS,  equa- 
tions (3.13),  RESIDUAL  FLEXIBILITY,  equations  (3.19),  and  MODAL  TRUNCA- 
TION, equations  (3.20),  formulations  each  consist  of  only  6+r  equations 
where  r is  selected  by  the  analyst.  The  value  of  r usually  ranges  from 
2 to  30.  The  simplest  formulation  of  the  equations  of  motion  is  the  QUASI 
STATIC,  equations  (3.4).  It  consists  of  only  6 equations  and  does  not 
explicitly  represent  the  elastic  deformations. 


The  development  costs  of  aircraft  systems  are  often  the  major  influ- 
ence on  the  selection  of  the  formulation  of  the  equations  of  motion  and, 
thus,  costs  must  be  weighed  in  any  analysis  as  an  additional  deciding 
factor.  In  the  event  that  two  or  more  of  the  six  formulations  for  the 
small  distrubance  motions  are  appropriate  to  a particular  flight  control 
analysis,  a cost-effectiveness  trade  can  be  constructed  similar  to  that 
presented  in  Table  4.2.  As  noted  in  the  table,  the  increase  in  accuracy 
requires  an  increase  in  analysis  cost  due  to  increasing  equation  complexity. 
The  reasons  for  this  increase  in  cost  and  complexity  are  that  both  com- 
puter program  run  time  and  program  size  increase  as  the  flight  control 
analysis  moves  from  the  QUASI  STATIC  to  the  EXACT  formulation.  The  pri- 
mary reason  for  the  increase  in  run  time  and  cost  is  the  necessity  to 
employ  highly  accurate,  computerized  routines  to  calculate  both  the  real 
and  complex  eigenvalues  and  eigenvectors  of  each  system  of  equations. 

In  addition,  the  RESIDUAL  FLEXIBILITY  formulation  requires  a non-standard 
flexibility  matrix,  the  "fret-free  flexibility  matrix",  that  is  not 
required  in  other  formulations.  This  flexibility  matrix,  related  to  the 
flexibility  matrix  of  the  cantilevered  aircraft,  must  be  recalculated, 
using  the  equation 

each  time  the  center  of  mass  is  changed,  C^>  , or  the  inertia  of  the  air- 
craft, M and  mj_,  is  altered.  These  two  changes  are  common  parametric 
studies  in  most  flight  control  analyses  and  are  particularly  important 
in  those  studies  that  are  used  to  design  and  analyze  the  reduced  static 
stability  augmentation  systems. 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


This  paper  has  related  the  common  formulations  of  the  equations 
describing  the  motions  of  elastic  aircraft.  The  conclusions  of  this 
study  are  that: 

1.  There  is  no  formulation  that  is  appropriate  to  all  flight  control 
analyses  of  elastic  aircraft.  The  large  disturbance  formulation,  equations 
(2.10),  must  be  carefully  specialized  for  each  analysis  due  to  their  im- 
mense complexity.  The  small  disturbance  equations  cf  motion  that  are 

most  appropriate  are  those  named  the  EXACT  formulation,  equations  (2.19). 
However,  in  the  interest  of  reducing  analysis  costs,  the  five  approximate 
formulations  of  section  3.0  may  also  be  appropriate.  In  no  case  should 
the  QUASI  STATIC,  equations  (3.4),  or  the  MODAL  TRUNCATION,  equations  (3.20), 
approximate  formulations  be  chosen  a priori. 

2.  The  selection  of  an  inappropriate  formulation  cf  the  equations 
of  motion  can  lead  to  large  errors  in  the  design  of  the  flight  control 
system.  The  actual  numerical  error  cannot  be  determined  at  this  time  due 
to  a lack  of  extensive  numerical  evaluations  of  the  terms  that  are  neg- 
lected to  develop  each  formulation.  The  numerical  evaluation  of  the 
significant  terms  neglected  in  each  of  the  formulations  should  be  a high 
priority  research  project.  Contemporary  aircraft  should  be  used  as  the 
study  configurations  in  this  research. 

3.  The  implementation  of  the  equations  of  motion  on  flight  simulators 
restricts  the  equations  to  the  RESIDUAL  STIFFNESS,  equations  (3.13),  the 
RESIDUAL  FLEXIBILITY,  equations  (3.19),  the  MODAL  TRUNCATION,  equations 
(3.20),  or  the  QUASI  STATIC,  equations  (3.4),  formulations.  This  restric- 
tion occurs  because  the  number  of  equations  must  be  small  such  that  the 
storage  capacity  cf  flight  simulators  (usually  less  than  40,000  central 
memory  locations)  is  not  exceeded.  Obviously,  the  QUASI  STATIC  formula- 
tion, consisting  of  only  6 equations,  is  the  simplest  to  use.  Usually 
this  formulation  does  not  use  all  of  the  central  memory  of  the  flight 
simulator  and,  thus,  the  remainder  of  the  central  memory  can  be  used  to 
represent  the  nonlinear  elements  of  the  control  system  or  the  nonlinear 
aerodynamic  effects  of  the  elastic  aircraft.  However,  the  QUASI  STATIC 
formulation  does  not  represent  the  dynamics  of  the  elastic  structure  and, 
for  these  analysis  cases,  the  other  three  formulations  must  be  applied. 

4.  The  formulations  that  may  most  easily  include  the  effects  of 
nonlinear  aerodynamics  are  the  RESIDUAL  STIFFNESS,  RESIDUAL  FLEXIBILITY', 

MODAL  TRUNCATION,  and  QUASI  STATIC  formulations.  This  is  because  these 
four  formulations  can  be  reduced  in  size  until  central  memory  storage  space 
in  the  digital  or  analog  computer  becomes  available  for  the  storage  of 

the  significant  nonlinear  aerodynamic  data.  The  penalty  for  the  reduction 
is,  of  course,  a reduction  in  accuracy. 
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